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ISOCAM DEEP SURVEYS UNVEILING STAR FORMATION IN THE MID-INFRARED 
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ABSTRACT 



Before having exhausted the helium in its tank on 
April 8, 1998, the Infrared Space Observatory had 
time to perform several complementary surveys at 
various depths and areas within selected regions of 
the sky. We present the results of some of the surveys 
done with the mid-infrared camera, ISOCAM, on- 
board ISO, and more specifically those surveys which 
were performed using the broad-band filter LW3 (12- 
18 /an). 

The counts obtained above 1 mJy by various surveys 
in the 12-18 jum band are in good agreement with 
the no-evolution model. On the contrary, a strong 
evolution is observed in the counts below 1 mJy by 
deeper surveys. This evolution seems to be due to a 
population of star-forming or post starburst galaxies, 
whose optical properties are not very different from 
the galaxies in the field. This implies that most of 
the star formation is hidden by dust and appears in 
the infrared as redshift rises. 



Keywords: cosmological surveys; mid-infrared; 
galaxy evolution. 



1. INTRODUCTION 



The so-called 'Madau plot' (Madau et al. 1996) 
which follows the rate of star formation per unit co- 
moving volume in the universe as a function of red- 
shift presents the advantage of being global enough 
to give an indication of the past history of star forma- 
tion as a whole, avoiding the difficulties encountered 
when trying to separate the respective role of indi- 
vidual galaxy types (Sp/Pec/E, hierarchical galaxy 
formation vs. individual collapse, duration of the 
bursts). However, based on optical-UV information 
only, it is now clear that this plot only sets a lower 
limit to the total star formation history of the uni- 
verse (see below). 

This plot also addresses the question of the heavy 
element content of the universe as a way to normal- 



ize the past history of star formation in the universe. 
Assuming that there was little margin left for more 
star formation than already found from optical-UV 
light one could see the lower limits given by Madau 
et al. 1996 as if they were firm constraints. Since 
then Lyman-break galaxies (also called UV drop-out 
galaxies) found by Steidel et al. (1996) appeared to 
be strongly affected by dust obscuration. Pettini 
et al. (1997) estimated that about two thirds of the 
emission due to star formation taking place in these 
galaxies is radiated in the infrared as indicated by 
the slope of their rest-frame UV spectral energy dis- 
tribution (see also Meurer et al. 1997, for the strong 
UV absorption at high redshifts). The metal content 
of the local universe presents enough uncertainties to 
leave margin for hidden star formation, especially at 
high redshifts where time flows faster per unit red- 
shift. Among those uncertainties, we shall quote the 
lack of galactic globular clusters with more than solar 
metallicity which would be required to calibrate the 
metallicity of metal-rich stellar populations found in 
elliptical galaxies; the presence of at least as much 
metals in the intra-cluster medium of galaxy clus- 
ters than in the galaxies themselves (Arnaud et al. 
1992, Mushotzky & Loewenstein 1997), with an over- 
abundance of a-elements which may indicate a biased 
IMF in the early stages of star formation in ellipti- 
cals (Elbaz et al. 1995); the still uncertain value and 
universality of the IMF. There is clearly space for a 
factor 2-4 in the heavy elements content of the local 
universe which implies a large freedom in the star 
formation history of the distant universe where time 
is so short. 

Although these indirect arguments already favour 
dust-enshrouded star formation with respect to 
optical-UV estimates in distant galaxies, the correc- 
tion for dust obscuration remains controversial and 
direct observations in the infrared are critical at this 
stage. The obvious way to go for it is the far infrared 
(FIR) since the spectral energy distribution (SED) 
of dusty star forming regions peaks around 60 //m. 
IRAS was not sensitive enough to probe the universe 
above z~0.2, except for some exceptionnally bright 
objects such as 10214+4724 (Rowan-Robinson et al. 
1991). FIRAS on-board COBE detected a diffuse far- 
infrared to sub-mm background (Puget et al. 1996, 
Guiderdoni et al. 1997, Fixsen et al. 1998). The ori- 
gin of this Cosmic InfraRed Background (CIRB) re- 
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mains controversial and could either be due to a local 
population of galaxies (z ~ 1) as proposed by Fall, 
Chariot & Pei (1996) or by a population of galax- 
ies distributed over a wider redshift range as sug- 
gested by Puget et al. (1996 and this conference, 
and Lagache et al. 1998). The recent detection of 
a population of distant dusty galaxies at 450 and 850 
/im with the sub-mm array camera SCUBA on the 
JCMT (Smail, Ivison & Blain 1997, Blain et al. 1998, 
Cimatti et al. 1998, Barger et al. 1998a, Hughes 
et al. 1998) tends to favour a shift of the peak of 
the star formation rate per unit comoving volume 
towards larger redshifts. These observations would 
favour the distant galaxy population option to ex- 
plain the CIRB. However, finding counterparts with 
spectroscopic redshift to these SCUBA galaxies will 
be the next step to answer this question. There will 
also still remain an uncertainty on the FIR part of the 
SED since sub-mm wavelengths only probe the right 
side of the peak, even at z=3 where the peak of the 
SED is rcdshiftcd to 240 /im. Reliable quantitative 
estimates of the star formation rate (hereafter SFR) 
will only become available when the rest-frame FIR 
part of the SED will be sampled. For this purpose, 
ISOPHOT on-board ISO, is not sensitive enough, al- 
though tantalizing results are already obtained at 175 
/im (Kawara et al. 1998, Lagache et al. 1998, Puget 
et al. 1998). The satellite FIRST will certainly play 
a major role in this respect. 

However, one can already see that hidden star forma- 
tion plays a major role. Although the mid-infrared 
range is less directly correlated with star formation 
than larger wavelengths, we will show that with a 
1000 times better sensitivity and 60 times better 
spatial resolution than IRAS, ISOCAM mid-infrared 
deep surveys confirm that most star formation is hid- 
den by dust and that this effect increases with red- 
shift (z < 1.5). 



2. ORIGIN OF THE MID-INFRARED 
EMISSION AND ITS LINK WITH STAR 
FORMATION 



The rest frame mid-IR emission of galaxies can be 
divided into three components (Puget & Leger 1989, 
Desert et al. 1990): 

• UIBs: the large Unidentified Infrared Bands 
(UIBs), detected at 6.2, 7.7, 8.6, 11.3 and 12.7 
/im as well as their underlying continuum, dom- 
inate the mid-IR emission below 12 /im (see fig- 
ure 1). The carriers of these UIBs contain Poly- 
cylic Aromatic Hydrocarbons (PAHs, Leger & 
Puget 1984, Puget & Leger 1989). The fact that 
in the wavelength range where the emission is 
dominated by UIBs, it globally scales with the 
intensity of the radiation field with only small 
changes in the spectrum shape is a strong indica- 
tion that the carriers of the UIBs are transiently 
heated by the absorption of individual photons 
(Boulanger 1998). The mid-infrared emission of 
galaxies does increase with the star formation 
rate though saturating for strong radiation fields 
(Boselli et al. 1998). 

• Warm dust (T>150 K): continuum at A > 
10 lira from Very Small Grains (VSGs) of dust 



(Desert et al. 1990). 

• Forbidden lines of ionized gas: Nell (12.8 /im), 
Nelll (15.6 fim), SIV (10.5 /im), Aril (7 /im). 
These lines are good indicators of the star forma- 
tion activity, particularly the Nelll/Nell ratio, 
when using the low resolution spectro-imagery 
mode of ISOCAM for nearby galaxies, but their 
contribution to large bands such as those used 
for the deep surveys is negligible. 



All three components affect the 12-18 /im ISOCAM 
filter (LW3) due to K-correction for galaxies between 
redshifts and 1.5. Indeed, the LW3 band includes 
emission from very small grains for low redshifts and 
becomes more and more contaminated by UIBs with 
increasing redshift. 
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Figure 1. SEDs from ISOCAM Circular Variable Fil- 
ter. Upper plot: Arp 220 (Charmandaris et al. 1998), 
UIBs nearly absent, high ratio of the 12-18 [im over the 
5-8.5 fim band. The displacement of the 12-18 \im (LW3) 
bandpass with redshift has also been plotted. Lower plot: 
M82 (Tran 1998), strong UIB features and low ratio of 
the 12-18 fim over the 5-8.5 fim band. 

The Far Infrared (FIR) flux is usually assumed to be 
correlated to the SFR once the cirrus emission is re- 
moved (see Scoville & Soifer 1991, Kennicutt 1998). 
The MIR emission is also a strong indicator of star 
formation, but the derivation of a star formation rate 
exclusively from the MIR flux is tricky: while many 
active star forming galaxies exhibit comparable in- 
frared SEDs, some are very different. The mid-IR 
over far-IR ratio can vary from less than 1 % to a 
few times 10 %. This results directly from the com- 
plex underlying physics of the mid-infrared emission 
which differs from the FIR one. For example, let us 
compare Arp 220 and M82 (see table 1 and figure 1): 

Both galaxies emit approximately the same 15 /im 
flux relative to visible light (respectively 3 , for Arp 
220, and 4, for M82, times more energy is radiated at 
15 /xm than at 0.55 /im). The contribution of the far 
infrared part of the spectrum relative to the MIR is 
substantially different: the ratio of the FIR to MIR 
light is 10 times larger for Arp 220 than for M82. This 
difference of an order of magnitude in FIR luminosity 
would imply one order of magnitude difference in star 
formation rate (using the formula given by Scoville & 
Soifer 1991 or Kennicutt 1998). Hence, the mid-IR 
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Table 1. Comparison of M82 & Arp 220: two tem- 
plate nearby starburst galaxies (H a —75 km.s -1 .Mpc -1 , 
distances corrected for the Virgo infall: d(M82)=3 Mpc, 
d(Arp220)= 73 Mpc). See Genzel et al. (1998) for a 
discussion showing that the emission of Arp 220 is dom- 
inated by its starburst activity and not its central AGN. 



The spectral energy distribution at the position of 
this spot presents two emission lines, Nell and Nelll, 
which do not appear as clearly on the brightest op- 
tical spots corresponding to the two galaxy nuclei, 
with a ratio of Nelll over Nell close to one, hence in- 
dicating the presence of massive young stars (figure 3, 
Mirabel et al. 1998). 
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flux cannot be considered as a direct measure of the 
star formation rate. However, Arp 220 is an extreme 
case and most galaxies behave more like M82 than 
Arp 220, so that the uncertainty is probably closer 
to 5 than to 10. 
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Figure 3. CVF spectra of 3 regions in the Antennae galaxy 
(NGC 4038/4039): the 2 bottom plots correspond to the 
2 galaxy nuclei, while the upper curve is associated to the 
brightest 15 fim spot in the overlapping region of the 2 
qalaxies. 
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Figure 2. mid-infrared image of the Antennae galaxies. 
The brightest spot lies in the overlapping region of the two 
interacting galaxies NGC 4038/4039, while the optical- 
UV emission peaks at the position of the two galactic nu- 
clei (Vigroux et al. 1996) 



3. ISOCAM DEEP SURVEYS 



Several extragalactic surveys have been performed 
using ISOCAM on-board ISO ranging from very large 
and shallow ones (ELAIS, P.I. M. Rowan- Robinson) 
to smaller and deeper ones (figure 4). 

As shown by the number counts plot (section 3.3.), all 
these surveys are complementary since they explore 
the behaviour of galaxies in the mid-IR over nearly 
4 orders of magnitude in the Log N-Log S plot. All 
surveys have been performed in regions of the sky 
relatively clean of cirrus and zodiacal emission: the 
Lockman Hole and Marano Field regions, but also 
the Hubble Deep Field region. The main limitation 
of these surveys comes from the presence of cosmic 
ray impacts on the detector for the detection and 
completeness limits (see below) and from the tran- 
sient behaviour due to the very low temperature of 
the detectors for the photometric accuracy. 



If the link between mid-IR fluxes and star forma- 
tion is not simple and implies the need for a multi- 
wavelength approach to reach reliable quantitative 
estimates of the star formation rate, the case of the 
Antennae galaxy is sufficiently explicit by itself to 
prove that hidden star formation can be unveiled 
through its mid-IR emission (Mirabel et al. 1998). 
ISOCAM mid-infrared imaging indeed revealed a 
bright spot (see figure 2, from Vigroux et al. 1996) 
which emits alone 15 % of the total mid-IR emission, 
within the overlapping region of the two interacting 
galaxies NGC 4038/4039 (Mirabel et al. 1998). This 
spot is clearly optically obscured and most of the op- 
tical emission arises from the two galactic nuclei. 



3.1. Data reduction 

ISOCAM data are subject to standard gaussian noise 
(photon and readout noises) and to errors associated 
with the flat-fielding and dark current substraction. 
But the main limitation of ISOCAM deep surveys 
comes from its thick and cold pixels: 

• because they are thick, ISOCAM pixel detectors 
are very sensitive to cosmic ray impacts (4.5 pix- 
els receive a glitch per second) . The behaviour of 
these glitches can be divided into three families: 



ISOCAM Extragalactic Surveys 




Figure 4- ISOCAM extragalactic surveys. 



— "normal glitches": the more common ones, 
which correspond to electrons and last only 
one or two readouts. They are easily re- 
moved with a median filtering (the combi- 
nation of several scales for the median al- 
lows the best correction). 

— "faders": these glitches as well as the follow- 
ing ones are probably associated with pro- 
tons and alpha particles. They induce pos- 
itive peaks in the detector response, which 
can last several readouts. Since ISOCAM is 
best used in the raster mode, a real source 
will look like these glitches, i.e. a positive 
response over the number of readouts spent 
on a given position of the sky. 

— "dippers": some glitches are followed by a 
trough extended over more than one hun- 
dred readouts. 

We have developed a tool based on multi- 
resolution wavelet transform which finds and re- 
moves the "faders" and "dippers". This tool, 
named PRETI, for Pattern REcognition Tech- 
nique for Isocam data (Starck et al. 1998), looks 
for the signal above the n-sigma level at differ- 
ent frequencies or scales in the wavelet trans- 
form of the original signal. Being extended, these 
glitches will be detected over several successive 
scales as typical patterns. Using their wavelet 
coefficients, one can substract a smooth function 
corresponding to the glitch in the raw signal and 
use the output signal in the co-addition of all pix- 
els of the raster that fall on the same sky position 
(see figure 5). 

• ISOCAM pixels are cold, so that electrons move 
very slowly within them and therefore induce 
a transient behaviour: a pixel will take several 
hundred readouts to stabilize when moving from 
the background to the position of a source on the 
sky and inversely. Because of time limitation, 
one is therefore limited to non stabilized signals 
which results in an uncertainty on the photome- 
try. This uncertainty is strongly reduced by the 



partial correction of the transient behaviour and 
by the use of simulations to define a statistical 
distribution of measured fluxes for any given in- 
put flux. The final uncertainty on ISOCAM deep 
survey fluxes is on the order of 20 %. 



A 



Figure 5. These plots show the decomposition of the sig- 
nal measured by a given pixel of the camera as a function 
of the number of readouts (Starck et al. 1998). From 
top to bottom: the first plot presents the typical time be- 
haviour of a given ISOCAM pixel. One can see from left 
to right: a "normal" glitch, a faint source which extends 
over readouts, a "dipper" , with long trough, and a 
"normal" glitch again, at the end of the trough. The sec- 
ond plot, shows the high frequency signal corresponding to 
the "normal" portion of the glitches. The third plot shows 
the reconstructed signal of the dipper's trough, after com- 
bining the information from several frequencies or scales 
of the signal. The fourth plot presents a gaussian at the 
position of the source with the extent of the source (we do 
not extract the signal at this frequency since it contains 
sources, but this plot is only intended to indicate the posi- 
tion of the source). The fifth plot shows the low frequency 
part of the signal, i. e. its baseline. The last plot presents 
the residual of the signal after extraction of the glitches 
and baseline. The source was substracted to show only the 
gaussian noise residual, due to photon and readout noise. 
Very faint sources are not visible on the time history of 
a given pixel and only appear after co-adding the residual 
of all pixels which pointed towards the same sky direction 
on a final raster map. 



In order to facilitate the separation of sources from 
cosmic ray impacts, ISOCAM surveys were per- 
formed using the raster mode with a redundancy 
(number of different pixels falling successively on a 
given sky position) ranging from 2 for the shallow- 
est survey (ELAIS) to 88 for the deepest surveys 
(Marano Field Ultra-Deep survey; 64 for the HDF 
field). 
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3.2. Simulations 

We performed simulations in order to quantify the 
sensitivity limit (minimum detected flux, below the 
completeness limit), the completeness limit (flux 
above which all sources are detected) and the photo- 
metric accuracy. These simulations were performed 
with real datasets (in order to simulate realistic 
glitches) in which we introduced fake sources includ- 
ing the PSF and their modeled transient: a long star- 
ing observation (more than 500 readouts) was used 
to simulate a raster observation. Real sources would 
extend over the whole observation and would be re- 
moved as a low frequency component of the signal 
while in a raster observation sources only remain dur- 
ing the number of readouts spent on a given position 
of the sky. Figure 6 shows as an example a sim- 
ulation of what an HDF-like ISOCAM deep survey 
would look like at 15 /im if the source distribution 
followed the model developed by Franceschini et al. 
(1991) with strong evolution (k — 3). The next sec- 
tion will present the real dataset to be compared with 
these simulations. 




Figure 6. simulated ISOCAM-HDF: from upper-left to 
bottom-right. The first image shows the distribution of 
sources expected from the model with strong evolution of 
Franceschini et al. 1994, where the number of galax- 
ies per unit luminosity per unit redshift is given by: 
N(L,z)=N(L,z—0)x expfnr) where k is the evolution pa- 
rameter (n=3)and T is the lookback time= (1q — t)/to. 
Fluxes are ranging from 0.1 /J,Jy to 1 mjy. Only the 
PSF is included and there is no noise associated with this 
"ideal" image. In the second image, the ISOCAM gaus- 
sian noise resulting from the readout and photon noise 
was added to the previous image. The last image presents 
the final image that ISOCAM would see including the 
noise introduced by cosmic ray impacts, shown in the third 
image. 

Finally, we also developed two independent data 
analysis techniques at Saclay (PRETI, Starck et al. 
1998) and Orsay (Tripple beam-switch technique, see 
Desert et al. 1998), in order to test the confidence 
level of the results. We find an agreement at the 20 % 
level corresponding to the photometric error of both 



techniques, an astrometric error close to the pixel size 
and compatible source lists. In the case of the Hub- 
ble Deep Field, we compared our results to the new 
results of the team of M. Rowan-Robinson (Rowan- 
Robinson et al. 1997, and references therein), which 
differ from their first analysis in the photometry at 
15 ^tm, because the techniques on each side were im- 
proved with time (see Aussel et al. 1998). 



3.3. RESULTS 



Source detection was performed on the output im- 
ages of PRETI (see figure 7 & figure 8 for the ISO- 
CAM HDF and Lockman Hole Deep Survey fields) 
using a spatial multi-scale wavelet transform. Hence 
multi-scale transforms were performed at two stages 
of the data analysis: first, on the temporal history of 
the detector pixels in order to remove glitches, then 
on the final co-added image for the source detection. 
The photometry of the Tripple Beam-Switch tech- 
nique from Desert et al. (1998) and PRETI are con- 
sistent at the 20 % level. 



ISO-HDF 15/j.m 




CENTER: R.A. 12 36 48.95 DEC +62 12 58.6 EOUINOX: J2000 

Figure 7. ISOCAM 15 fim HDF image from Aussel et al. 
(1998). Compare this image to the simulation shown in 
the previous image. The excess of detections found in the 
real image reflects the high evolution paramater needed to 
reproduce the number counts below (k= 6 instead of 3). 

Several hundreds of simulated images were done in 
order to estimate from a Monte Carlo statistical ap- 
proach the limits of our results in terms of photomet- 
ric accuracy and completeness. We have included this 
information in the Log N-Log S curve (figure 9): error 
bars on the fluxes and the flux range where complete- 
ness is achieved result directly from the Monte Carlo 
simulations. 

The number counts are presented here in the two 
classical forms: the log N-log S (figure 9), gives the 
number of objects which flux is higher than a given 
flux limit S(mJy), and the differential counts nor- 
malized to euclidean (figure 10). No evolution would 
give a slope of 3/2 in the Log N-log S in a euclidean 



6 



Lockman Hole Deep Survey 12-18 microns (LW3) 
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Figure 8. ISOCAM 15 pm Lockman Hole Deep Survey 
image from Elbaz et al. in prep. 



space (2.5 in the differential counts, by derivation). 
In a relativistic universe the slope of the no evolution 
curve is even lower than that as shown by the no- 
evolution curve from Franceschini (1998) overplotted 
on figure 9. 
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Figure 10. euclidean normalized differential counts: 
ISOCAM-HDF (diamonds), Lockman Hole deep survey 
(stars), ELAIS (crosses). The models are the same as in 
the log N-log S plot. 



-10 12 
log ]0 (S (mJy)) 



ISOCAM Deep Surveys 12-18 pm 



-6 



Hubble Deep Field - slope: 1.3 

ckman Hole Deep Survey - slope: 2.1 

Lockman Hole Shallow Survey" 
ELAIS - slope: 1.3" 



Models (Franceschini): 

no evolution 

evolution 

new population 



-2 



-1 1 

log, S (mJy) 



Figure 9. log N( S)- Log S: number of objects, N, detected 
above a given flux level, S(mJy), with ISOCAM in the 
15 pm broad band filter LW3. The data are consistent 
although they come from 5 different origins. From high 
fluxes to low fluxes one finds: ISOCAM-HDF (from Aus- 
sel et al. 1998), Lockman Hole deep survey (Elbaz et al, 
in prep), Lockman Hole shallow survey (Elbaz et al. in 
prep), ELAIS (preliminaty results from the ELAIS con- 
sortium, P.I. M. Rowan-Robinson, Oliver et al. 1998), 
IRAS 12 pm (Spinoglio et al. 1995). Dashed line: no 
evolution, normalized to IRAS and including PAH fea- 
tures. Dashed-dotted line: strongly evolving population in 
N(L,z)=N(L,z=0)x exp(nr) with k=6. Plain line: sum 
of the no evolution + strongly evolving population. The 
strongly evolving population becomes dominant below 1 
mJy 



Figure 11. euclidean normalized differential counts: same 
as above. Models of Xu et al. (1998) normalized to the 
same IRAS point as in figure 10 . Plain line: pure lu- 
minosity evolution, p oc (1 + z) A . Dashed-dotted line: 
pure density evolution, p oc (1 + z) 3 . Dashed line: no 
evolution. 



3.4. Interpretation 

For the first time, we present number counts over 

4 orders of magnitude ranging from 100 /iJy, in the 
ISOCAM-HDF (Aussel et al. 1998), to 50 mJy, in 
the ELAIS fields (preliminary results from the ELAIS 
consortium, P.I. M. Rowan-Robinson, Oliver et al. 
1998) and up to 300 mJy, including the IRAS 12 
pm counts. It can be seen on both the log N-log 

5 (figure 9) and the differential counts (figure 10) 
that above the 1 mJy level, the counts follow the no- 
evolution model, while below this flux level, a rapid 
rise appears (clearly seen on the differential counts 
plot). Such a trend cannot be reproduced either by 
the PAHs features, which effect on the K-correction 
is too weak, or by a continuous evolution such as : 



1. pure luminosity evolution as described in 
Franceschini et al. (1991): N(L,z)=N(L,0) x 
exp(rcr) where k is the evolution parameter and 
r = (to—t)/to is the lookback time. No value of k 
allows to fit the data, except if we consider a sep- 
arate population of objects which becomes dom- 
inant below the 1 mJy level, i.e. above z~0.4. 

2. pure luminosity evolution used by Pearson & 
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Rowan-Robinson (1996) to fit the IRAS 60 fan 
number counts, proportional to (1 + z) 3 and cal- 
culated by Xu et al. (1998) at 15 (im (figure 11). 

3. the pure density evolution, in (1 + z) 4 , from Xu 
et al. (1998, figure 11). 



• a few bright infrared objects make most of the 
cosmic background, since the few detections at 
15 fan. in the HDF radiate more energy ,when 
integrated over the whole infrared range, than 
the large number of optically detected galaxies. 
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Figure 12. Histogram of the I-K rest-frame colours 
of all HDF+ Flanking fields galaxies (dashed line) and 
HDF+ Flanking fields galaxies detected at 15 fim above 
100 flJy. The rest-frame colours were calculated using 
the HDF and Flanking Field catalog of Barger et al. 
(1998b) with a method similar to the one used by Lily 
et al. (1995). The I band K-correction was computed by 
us using the Coleman et al. 1980 SED and the K band 
K-correction is from Cowie et al. (1994). The difference 
in the median colours is only 0.2 dex. 

Extreme evolution needs to be added to the no- 
evolution curve, which could be associated with a 
galaxy population dominating the counts below 1 
mJy, i.e. above z=0A (in this model). The evo- 
lution parameter, k, required for this "population" 
is k=6, while it was 3, in the original model with 
strong evolution of Franceschini et al. (1991). This 
rapid increase of the number counts below 1 mJy, 
corresponds to a high flux density in the mid-IR 
going towards faint fluxes. Integrating the num- 
ber counts over the whole flux range one finds that 
about 2 x 10 -9 W.mT 2 .sr _1 is emitted at 15 fan and 
3 x 10 -9 W.m~~ 2 .sr" 1 when extrapolating the counts 
assuming no further increase in the slope. This cor- 
responds to 30-45 % of the cosmic background seen 
with the HST in the I-band (7 x 10~ 9 W.mT 2 .sr" 1 ). 
Locally about one third of the energy is radiated in 
the infrared (Soifer & Neugebauer 1991) and mainly 
in the FIR since the FIR over MIR ratio is typically 
3 for normal galaxies and stronger than 10 for star- 
bursting galaxies. The cosmic infrared background 
is therefore larger than the optical background when 
integrating over a larger redshift range than IRAS, 
which implies that going towards higher redshifts: 

• dust obscuration is stronger at higher redshift. 
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Figure 13. Fraction of Peculiar versus Spiral galaxies in 
the HDF+ Flanking fields for galaxies detected at 15 fim 
(plain lines) and all galaxies (dashed lines). Above Z—0.4, 
peculiar galaxies dominate the distribution of galaxies de- 
tected with ISOCAM while optically selected galaxies are 
mainly spirals. This population of peculiar galaxies may 
produce the rapid evolution below the 1 mJy level found 
in the 15 fim number counts. 

This result agrees with SCUBA data (Barger et al. 
1998a, Hughes et al. 1998) and with the CIRB es- 
timates (see Puget et al., these proceedings). Ad- 
ditional conclusions can be reached only after iden- 
tification of the sources, and multi- wavelength and 
redshift measurements. We are currently gathering 
such data for the sources in the Lockman Hole. In 
the HDF and in one of the CFRS fields, surveyed 
with ISOCAM by Flores et al. (1998), such data are 
available. We summarize some of the main conclu- 
sions obtained from this source by source comparison 
over ~ 60 sources within the HDF+flanking fields 
with redshifts in the range 0.4-1.3, as well as in the 
CFRS : 

• at z>0.5, peculiar, mergers and post-starburst 
galaxies are much more frequent in the ISO base 
than in the optical base (Flores et al. 1998, Elbaz 
et al. (in prep.), see figure 13). It is likely that 
the additional population required to account for 
the number counts is made of such galaxies. 

• Flores et al. (1998) attempt to fit the multi- 
wavelength data over 10 sources for which radio 
data are available. In this way, they estimate 
that the universal star formation at z~l is en- 
hanced by a factor 2.3 with respect to previously 
estimated UV fluxes. 

Using available data on the HDF, especially optical 
and NIR data from Barger et al. (1998b), we find 
that galaxies detected above 100 /zJy (completeness 
limit reached in the HDF, Aussel et al. 1998) do not 
show any clear signature in their optical colours (see 
figure 12, Elbaz et al. in prep). The I-K colour of ISO- 
CAM detected galaxies is only 0.2 dex redder than 
all HDF galaxies. 
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3.5. CONCLUSION 



We showed that even with very sensitive telescopes 
like the HST, optical observations need to be com- 
bined with infrared imaging in order to select the 
major events of star formation in the universe. The 
sensitivity required for the NGST to detect all HDF 
optical galaxies, i.e. 2xl0 6 galaxies per square seg- 
ree, at 15 (im would be about 0.02 fi3y if the slope 
flattens as shown by the model on figure 9 and 0.3 
/xJy if a moderate evolution with a slope close to 
1.3 continues below the 100 ^Jy level. For a large 
fraction of these galaxies, the NGST would resolve 
the regions of star formation in actively star form- 
ing galaxies and change our view on galaxy evolution 
without suffering from dust obscuration. 



REFERENCES 



Arnaud M., Rothcnflug R., Boulade O., Vigroux L., 
Vangioni-Flam E. 1992, A&A 254, 49 

Aussel H., Cesarsky C, Elbaz D., Starck, J.L. 1998, 
to appear in A&A 

Bargcr A. J., Cowie L. L., Sanders D. B., Fulton E., 
Taniguchi Y., 1998a, Nature 394, 248 

Barger A. J., Cowie L. L., Trentham N., Fulton E., 
Hu E. M. et al., 1998b, AJ submitted, available at 
http://wwvj. if a. hawaii. edu/ cowie/hdflank/hdflank. html 

Blain A. W., Smail I., Ivison R. J., Kneib J.-P. 1998, 
submitted to MNRAS, astro-ph/9806062 

Bosclli A., Lequeux J., Sauvage M., Boulade O., 
Boulanger F. et al. 1998, A&A 335, 53 

Boulangcr, F. 1998, Star Formation with the Infrared 
Space Observatory, Lisbon, Portugal, edited by Joao 
Yun and Rene Liseau, Publisher: ASP, as volume 
132 of ASP Conference Series, p. 15. 

Charmandaris, V. et al. 1998, in preparation 

Cimatti A., Andreani P., Rottgering H., Tilanus R. 
1998, Nature 392, 895 

Coleman G. D., Wu C. C, Weedman D. W. 1980, 
ApJ Suppl.Ser. 43, 393 

Cowie L. L., Gardner J. P., Hu E. M., Songaila A., 
Hodapp K.-W., et 1994, ApJ 434, 114 

Desert F.-X., Boulanger F., Puget J.-L. 1990, A&A 
237, 215 

Desert, F.-X., et al. 1998, to appear in A&A 

Elbaz, D., Arnaud, M., Vangioni-Flam, E. 1995, 
A&A 303, 345 

Elbaz, D., et al. 1998, in preparation 

Fall S.M., Chariot S., Pei Y.C. 1996, ApJL 464, 43 

Fixsen D. J., Dwek E., Mather J. C, Bennett C. L., 
Shafcr R. A. 1998, ApJ, in press 

Flores H., Hammer F., Desert F.-X., Cesarsky C. J., 
Crampton D. et al. 1998, to be submitted to ApJ 

Franceschini A., Toffolatti L., Mazzei P., Dancsc L., 
De Zotti G., 1991, A&A Suppl.Ser. 89, 285 

Franceschini, A. 1998, in T.Thuan, G. Mamon, and 
J.T.T.Van (eds), Extragalactic Astronomy in the In- 
frared, Proceedings of the XVIIth Moriond Astro- 
physics Meedting, Ed. Frontiere, p. 509. 



Genzel R., Lutz D., Sturm E., Egami E., Kunze D. 
et al. 1998, ApJ 498, 579 

Guidcrdoni B., Bouchet F.R., Puget J.-L., Lagache 
G., Hivon E. 1997, Nature 390, 257 

Hughes D. H., Serjeant S., Dunlop J., Rowan- 
Robinson M., Blain A. et al. 1998, Nature 394, 
241 

Kawara, K., Sato Y., Matsuhara H., Taniguchi Y., 
Okuda H. et al. 1998, A&A 336, L9 

Kcnnicutt, R.C. 1998, ApJ 498, 541 

Lagache G. 1998, to be submitted 

Leger A., Puget J.L. 1984, A&A 137, L5 

Lilly S. J., Tresse L., Hammer F., Crampton D., 
Lefevre O., 1995, ApJ 455, 96 

Madau P., Ferguson H. C, Dickinson M. E., Gi- 
avalisco M., Steidel C. C, Fruchter A. 1996, MN- 
RAS 283, 1388 

Meurer G. R., Hcckman T. M., Lehnert M. D., Lei- 
therer C, Lowenthal J. 1997, AJ 114, 54 

Mirabel I.F., Vigroux L., Charmandaris V., Sauvage 
M., Gallais P. et al. 1998, A&A 333, LI 

Moorwod A. F. M., Lutz D., Oliva E., Marconi A., 
Nctzer H. et al. 1996, A&A 315, L109 

Mushotzky R. F., Loewenstein M. 1997, ApJ 481, 
L63 

Oliver S. J. et al. , 1998, IAP workshop on Wide- 
Field Surveys. 

Pearson C, Rowan-Robinson M., 1996, MNRAS 283, 
174 

Pettini M., Steidel C. C, Adelberger K. L., Kellogg 
M., Dickinson M. et al. astro-ph/9708117 

Puget J.-L. Leger, A. 1989, Ann.Rev.Ast.Ap. 27, 37 

Puget J.-L. , Abergel A., Bernard J.-P., Boulangcr 
F., Burton W. B. et al. 1996, A&A 308, L5 

Puget, J.-L. et al. 1998, these proceedings 

Scoville, N.Z., Soifer B.T. 1991, in Massive Stars in 
Starbursts, ed. Leitherer C, Walborn N., Heckman 
T., Norman C, Cambridge Univ. Press, p. 233 

Soifer B.T., Neugebauer G. 1991, AJ 101, 354 

Spinoglio L., Malkan M.A., Rush B., Carrasco 
L.,Recillas-Cruz E. 1995, ApJ 453, 616 

Starck J.-L., Aussel H., Elbaz D., Cesarsky C. J., 
1998, in T.Thuan, G. Mamon, and J.T.T.Van 
(eds), Extragalactic Astronomy in the Infrared, Pro- 
ceedings of the XVIIth Moriond Astrophysics 
Meedting, Ed. Frontiere, p. 509. 

Steidel C.C., Giavalisco M., Pettini M., Dickinson 
M., Adelberger K.L. 1996, ApJ 462, L17 

Rowan-Robinson M. , Broadhurst T., Oliver S. J., 
Taylor A. N., Lawrence A. et al. 1991, Nature 
351, 719 

Rowan-Robinson M., Mann R. G., Oliver S. J., Ef- 
stathiou a., Eaton N. et al, 1997, MNRAS 289, 
490 

Smail I., Ivison R. J., Blain A. W. 1997, ApJL 490, 
5 

Tran D. 1998, PHD thesis, Universite de Paris XI. 

Vigroux L., Mirabel I. F., Altieri B., Boulangcr F., 
Cesarsky C. J. et al. 1996, A&A 315, L93 

Xu C, Hacking P. B., Fang F., Shupc D. L., Lonsdale 
C. J. et al., 1998, ApJ in press, astro-ph/9806194 



